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DESIGN OF A NONLINEAR AOAPTIVE FILTER FOR 
SUPPRESSION OF SHUTTLE PILOT-INDUCEO 
OSCILLATION TENDENCIES 

John W. Smith and John W. Edwards 
Dryden Flight Research Center 


INTRODUCTION 


On the final flight of the space shuttle's approach and landing 
tests, a pilot- induced oscillation (PIO) was experienced during the flare 
and landing. The frequency of the PIO in the pitch axis was approx- 
imately 3.5 radians per second, and the pilot utilized a significant portion 
of his full control authority. Subsequent analysis confirmed the exist- 
ence of PIO tendencies in this configuration, and indicated that the 
source of the problem was a combination of basic shuttle handling 
qualities, time delay through the digital flight control system computers, 
and rate limiting of the elevator actuators. In addition, the cockpit of 
the shuttle vehicle is at the approximate center of rotation of the 
vehicle for pitch maneuvers, which may deprive the pilot of some motion 
cue information. 

High pilot gain during a critical task, coupled with the ability to 
rate limit surface actuators with large pilot inputs, is a significant 
factor in such PIO problems (refs. 1 to 3). Linear compensation tech- 
niques are of limited usefulness in solving this problem because of the 
phase lags they introduce. This report describes a PIO suppression 
(PICS) filter that was designed to reduce pilot gain in PIO situations 
with little or no phase lag. The filter is implemented as a nonlinear, 
adaptive element of the pilot's stick gearing schedule and operates by 
sensing the frequency and amplitude of the pilot's stick inputs. 

The primary effect of attenuating the command path gain is to reduce 
the amount of rate limiting. Reducing the gain near the crossover 
frequency without phase lag may also result in some improvement in 
handling qualities, as shown in reference 4. 


SYMBOLS AND ABBREVIATIONS 


A, B, C, D, 
E, X, Y 

ADEP 

AMP 


dummy variables 

2 

output of the smoothing filter, deg 
ampl itude 

normal acceleration, g 



BN output of the control activity estimator, deg/sec 

CPS samples per second 

DCGN output of the deadband nonlinearity, deg 

DEC suppressed output of the pitch stick shaping function, deg 

DEC 1 unsuppressed output of the pitch stick shaping function, 

deg 

DEP pilot input, deg 

EDEP ratio of filter parameters, per sec 

F(s) transfer function 

K c forward loop gain, per sec 

K(q) gain constant 

K pilot gain 

p 0 

M, N indexes 

NPTS number of points 

PiLN initial controller input equal to DEPN, deg 

PIO pi lot- induced oscillation 

PIOS PIO suppression 

RKP control activity estimation equation 

SKQ output of the F-2 schedule 

SQADEP square root of the filter output along the amplitude 

path, deg 

SQXPN square root of the filter output along the rate path, 

deg/sec 

s Laplace transform variable, rad/sec 

t (T) time, sec 

XPN squared and filtered output along the rate path, 

2 2 
deg /sec 
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pilot model 


2 2 

output of th. control activity estimator squared, deg /sec 
elevator deflection, deg 

pitch attitude, deg 
frequency, rad/sec 
phase angle, deg 

Subscripts: 

A/C aircraft 

e error 

A dot over a quantity denotes the first derivative with respect to time. 

SHUTTLE PITCH RATE COMMAND SYSTEM 


V s) 

ZN 

6 e 

e 

U) 

♦ 


A pitch rate conmand systea is used on the space shuttle for the 
approach and landing phases of flight. The systea is impleaented as a 
digital flight control systea resident in the shuttle's flight control 
coaputers. Figure 1 shows a conceptual block diagraa descriptive of the 
eleaents important to closed-loop pitch attitude control. The pilot, des- 
ignated as Y (s), applies a torque to a rotational hand controller which 
p 6 

coaaands an output given by the following pitch stick shaping function: 

DEC = (0.36 + 0.0484 |DCGN| )DCGN 


The input to the stick shaping fur -♦ion, DCGN, results froa operating on 
the pilot's stick position, DEP, with a ±1.15° deadband. The pitch rate 
command is K £ * DEC in (deg/sec )/deg; and the gain constant for the 

approach and landing, K 0 remains constant at a value of 0.4. Pitch rate 

control is achieved by positive feedback around the power actuator and 
servo, yielding an equivalent system for the actuator loop, with a "free s" 
in the denominator. As such, the error signal, which is made up of 
the commanded response and aircraft pitch rate, is continually driven to 
zero. The integration rate is proportional to the scheduled gain due to 
dynamic pressure, K(q), and is approximately equal to 1.446. The 
remaining terms left in the closed-loop transfer functions are: 


F(s) = 1.395 


(s/1.5 ♦ l)(s/l. 8 + 1) 

(s/21.5 t flKs7T7 ♦ 1) 
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Inside the actuator and servo loop, the rate at which the actuator can 
move is restricted to ±20 deg/sec. The filters, gains, summations, rate 
restrictions, and pitch stick shaping function are implemented in the 
shuttle's digital flight control computers. 


PILOT-INDUCED OSCILLATION SUPPRESSION FILTER 


Description and Design Logic 

In general, the objectives of the PIOS filter are simple: to modify 
the stick shaping function, DEC = f(DCGN), as a function of the frequency 
and amplitude of the pilot's input to acquire the desired gain reduction; 
and, at the same time, to minimize phase lags introduced by the filter. 

The filter developed to achieve these objectives is diagramed in 
figure 2. The blocks of the diagram present all the elements, 
mathematical computations, and logic of the filter in a sequential manner. 
The dynamic elements, shown as Laplace transformed functions, are imple- 
mented as digital filters. The FORTRAN program used to implement the 
filter is listed in appendix A. 


The PIOS filter contains a frequency- and amplitude- sensitive 
element composed of a second-order lag-lead filter 

(i - ^a ); two rectifiers; two first-order filter smoothers; and 

4 (s + ior 

a modifying function, SKQ, for the quadratic term of the Stic shaping 
function. The PIOS filter operates by obtaining an estimate, EDEP, of 
the dominant frequency of the pilot's input. This estimate is used to 
modify the quadratic term of t. a stick shaping function in order to 
preclude the possibility of rate limiting the control surface actuators. 


The pilot's input, DEP, is displaced by the deadband (±1.15°) and 
then goes to both the pitch stick shaping function and the elements of 
the PIOS filter. The pitch stick shaping function implemented on the 
shuttle (fig. 3) is give'' by 

DEC 1 = (0.36 + 0. 0484 |DCGN| )DCGN 


The combination of linear and quadratic terms in the gearing is 
intended to give acceptable handling qualities for small inputs while 
also allowing full control surface authority. 
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The output of the summation junction in the PIOS filter paths, BN 
(fig. 2 and below), essentially serves as a control activity estimator. 

It is defined as follows: 

BN = [l - l OCGN 

L 4 <s ♦ 10) Z J 

= &(» ♦ 13.33 ) 

1 4(s + ior 

2 

This quantity, BN, then is rectified (X ) and smoothed by a low pass 
filter. The square root of the resulting quantity, SQXPN, is proportional 
to the frequency and amplitude of the low frequency components of the 
pilot's input. The sane mathematical operations are carried out along 
the amplitude-estimating path, providing a quantity proportional to the 
amplitude of the pilot's input. The F-l schedule (fig. 4(a)) 
serves two purposes. First and foremost, the schedule is necessary to 
avoid subsequent division by zero. However, to lessen the suppressed 
response for step inputs of 5° and less, the minimum value is limited to 
4.5°. A ratio of the two paths is obtained to eliminate the amplitude 
effect, making EDEP primarily proportional to the frequency of the control 
activity. The F-2 schedule (fig. 4(b)) controls the amount of suppression 
provided by the filter. The minimum value of the F-2 schedule (-0.15) 
limits the maximum suppression of the filter. The design philosophy 
regarding the time constants and schedules for this particular example 
is given in appendix B. 


Digital Computations and Simulations 

Figure 5 presents the results of a digital simulation using the 
FORTRAN listing givpn in appendix A. Pertinent filter elements are shown 
as a function of time when forced by a sinusoidal pilot input, DEP. The 
peak input amplitude is 10° at a frequency of 3 radians per second for 
3.25 cycles. The "notching" of BN and ZN is due to the deadband. However, 
this effect is smoothed out by the first-order filter. Parameters down- 
stream of the first-order smoothing filters illustrate the 3.3 second 
settling time. In particular, the SKQ trace is approaching its minimum 
value of -0.15. Fo r ‘ this set of conditions, the F-2 schedule was 
designed to almost limit the maximum amount of suppression and the final 
value of the stick shaping function, DEC. DEC is reduced to about 43 
percent of the unsuppressed output, DEC 1, and this value is almost 
reached after three cycles. The steady-sta*® suppression and settling 
time for this example were chosen as the design objectives (appendix B). 

Cross plots of the output-to- input relat \>nships for input amplitudes 
of 17°, 10°, and 5° are presented in figure 6. The sinusoidal input 
frequency is 3 •'dians per second and oscillates for five complete cycles. 
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For the 17° input (fig. 6(a)) the final suppressed output is approximately 
5°, and was nearly reached after two cycles. This particular character- 
istic of the PIOS filter is due to the time constant of the first- 
order filter and the F -2 schedule. The final suppressed output is about 
28 percent of the unsuppressed output shown in figure 3. There is very 
little hysteresis, even at the extreme input, and none at 0° of DEP. 

Figure 6(b) presents the cross plot for a 10° input. The same general 
trends as for the 17° input amplitude response are evident. The final 
suppressed output for a 10° input is 3°, which is about 43 percent of the 
unsuppressed output. Figure 6(c) shows a cross plot of the output-to- 
input ratio for 5° of input amplitude. Very little suppression is evident. 

The final output is about 75 percent of the unsuppressed value. This 
characteristic is due primarily to the minimum value (4.5) set on schedule F-l. 

Information similar to that presented in figure 6 can be summarized 
over the frequency range of interest. Figure 7 presents the final sup- 
pressed output of the PIOS filter for three different input amplitudes 
(17°, 10°, and 5°) over a frequency range from 0 to 7 radians per second. 

Two characteristics are apparent. First, at frequencies above 4 radians 
per second the suppressed output is constant for all conditions. This is 
because the minimum value in schedule F-2 (-0.15) has been reached. 

Second, at frequencies below 3 radians per second, the reduction in 
output is almost linear with frequency. Finally, at low input 
frequencies, ti-ere is very little suppression for an input amplitude of 
5°. 


Figure 8 presents the time response of the pilot command with and 
without the PIOS filter for a 5° and 10° step input. The trace for the 
5° step input shows the effects of the design consideration placed on the 
F-l schedule in that SQADEP had a minimum value of 4.5 (fig. 4(a)). The 
output of the PIOS filter, DEC f~r a 5° step is within 85 percent of the 
unsuppressed value. For a 10° :U! inout, the output of the shaping 
function varies from a ninimuro value oV 3.7° at 0.25 second to almost 6° 
at 4 seconds. The rise rate of the output response, DEC, is governed by 
the time constant (3.3 sec) in the first-order smoothing filter. This 
time-dependent response would reduce the amount of surface rate limiting 
due to large step inputs. 


Open- Loop Response 

The actual pilot input, DEP, during the shuttle landing when the PI0 
occurred was used as the input to the filter to illustrate the PIOS filter's 
behavior. Even though this simulated approach is strictly open loop, the 
output of the PIOS filter is indicative of the control surface motion to 
be expected under similar circumstances with the PIOS filter in operation. 
Figure 9 presents the response of pertinent PIOS filter parameters due to 
the pilot input, DEP, as generated from the data recorded during the 
flight. During the approac. from 0 to approximately 20 seconds, the 
amount of suppression varies, as is evident in the governing parameter 
SKQ. For this time span, SKQ is always less than 1.0 and averages about 
0.5. During the PI0, between 25 and 35 seconds, the peak inputs were 
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about ±10° at a frequency of between 3 and 4 radians per second. The SKQ 
parameter for this period of tine would have reached a mini nun value at 
about 25 seconds. This would have occurred just prior to landing and 
throughout the landing phase. The suppressed output, DEC, compared with 
the unsuppressed output, OEC 1, is noticeably less throughout the latter 
part of the flare and landing nanuever. For this example and with respect 
tc the predicted rate limiting, the PIOS filter is designed to offer the 
desired amount of suppression following the logic developed in appendix B. 


Closed- Loop Analysis 

Observation of the actual closed-loop behavior of the P70S filter 
would undoubtedly permit a more meaningful assessment of its worth. 
Therefore, a constant gain pilot model K (fig. 1) was mechanized on a 

p 0 

simulator of the shuttle flight dynamics. Simulator (fig. 1) responses 
were obtained with and without the PIOS filter included. Figure 10(a) 
shows the characteristics of the simulated shuttlecraft with pitch attitude 
feedback and basic stick shaping for a pilot gain. K , of 7.2. This gain 

p 0 

was selected for this example as the gain where instability would begin to 
occur. With the PIOS filter (fig. 10(b)), the system was stable, with a 
subsidence ratio of -0.38, and would have taken about one cycle to damp out. 
Both systems were pulsed with a pitch input of 1°. 

Figures 11(a) and 11(b) present the variation of subsidence ratio 
and frequency, respectively, with gain. At a pilot gain of 4, both systems 
are well damped. With the PIOS filter, about twice the gain margin is 
realized, and for higher gains the oscillations are bounded. The 
oscillation frequency approaches 3.0 radians per second. The damped 
frequency is the same for low pilot gain and somewhat less with the PIOS 
filter below the divergent gain. 


CONCLUSIONS 


A nonlinear adaptive filter was designed to suppress the pilot- 
induced oscillations (PI0) tendencies of the space shuttle vehicle. From 
digital computations and simulations, it can be concluded that: 

(1) The PIOS filter attenuates the gain of the pitch stick shaping 
function without adding phase lag to the system. 

(2) With the pitch attitude looo on a shuttle simulator model 
closed, the PIOS filter demonstrated a gain margin increase by a factor 
of about two. 


Dry den Flight Research Center 
National Aeronautics and Space Administration 
Edwards , California, March 5 , WHO 
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APPENDIX A - COMPUTER PROGRAM FOR PI LOT- INDUCED 
OSCILLATION SUPPRESSION WITH SINE WAVE INPUT 


P#OGRA* SFILTE* <!NPUT#OOTPLT#1AFEWNPU!,TAPF2) 

MARCL1ST /klLNAh/ AHP #OHEGA»CPS 

C ***** ifcT INITIAL AhO CONSTANT VALUES* 

CALL PLOT? (fc.-WJI 
CALL FAC TON (0*3937) 

i reap aaoo) 

IOC FORMAT (#10) 

IF UcFm.NE.OI G3 TO 99V 
BACKSPACE 1 
OF AO FILNA* 

PAINT 1*1# ANP#0*EGA#CPS 

1C* FOAHAT (1m x ,-ahp - •, F 9*1# IX, •PNECA - "# F3*1#5X# **CPS • "*F5.l*/) 
POINT 102 

*v2 FOft“AT(lM,A* f "TinE”»9X#”0tt'”#9X#-ftN*f9X#"ZN"#9X»”XPN"#9X»”S0XPN"» 
9 9l# •*SQADFP"#9X#“EDEP ,, #9 a# *SK9"#9X# "0§C1*#9X# "DEC"# /) 

T • 1.0/CPS 

NPTS • 9* 0*6*20 32*CPS /PNEGAaO*; 

A - — 2.0*FXP f — 2C*0*T I 
9 • tXP(-A0,0*T ) 

C ■ -2.o*tXP(-10*0*T) 

P ■ FXP|-23.**TI 
E-EXP(-*30*T> 

F«*3*U*-F) 

RKP • ( tl.O— fcXP (— 10*d*Tl )**2 )/ ( 1 1*0— EXP I— 2v*0*T ) )••?) 

VNN1 • TNN2 * OEPNKl • DEPNM2 - XPNN1 • ZNNl • 0*0 
A0EPN«At)EPNl*0fcPK2»CcPN21*b* 

C **•* LOJP TO CONFUTE LEP AND DEC FOR 5 CYCLES* 

N - -I 

> * • N*1 

IF (N. 6T.Nr TS ) GG TO A 
PILN ■ ANP ♦SlN(h*OHCGA*T ) 

OFPN*** ILN 

IF (ASS(utPN).lt. 1.15) TCGN *0*0 
IF I C OE?N ) *GT • J.1SI OCGN -DkPh-1.15 
IF UucPM.LT* -a* 1* I OCGN • OtPN ♦x.l? 

OtFK "OCGN 

YN • -w*YKH]-O*YNN2*OKP*(0EPN«A*UEPNN«*B*DFPNN2) 

PH* ( Dt P 4-»YN) 

m>tau>**? 

XPh«t**PNNl*FM IN*?**! ) 

S3XPN-S9PT(>Ph) 

OcPN # «C c PN*DFPN 

ADFPN«l*AOcFWl*f *(D t PN TEPN21) 

SO AC f P • (SlPTUOFFN) I 

IF (SOAD'P.LT.A.SI SOAOEP«A*3 

fcOFP-SOXPK/SCAOEP 

T*W»l.-<l./*<A>*fD e P 

IF (SRG.U.-.15) SKO --.15 

Oifl • r )C&NA( - .36*O.CAFA**oS(D r GN)) 

JFC-OCGN* C0.3t *SRO*C.CACA* AESUCGM > 

IF(VPT5.GT,1 : v *. c .ANO. rCO!N#I''U).NF.'U go TO 3 
IFINPT^CT.K** .ASC* "OCtN# a*M *Mf . M cr T n 3 
TINF • N* 

°P IN T 1^3#1 ilk . PiLN#BN#ZN# *PN»SQ*PN # SQACk P# t DFF,S^O,OcCl#rEC 
*w 3 FCOflAT( *X,U( ► W . 5 , 2 * ) ) 

* *RiTl (21 PlLNfUrc 

YSN 2 - f N *1 
TN *1 - YN 
0EMW2 • oepNni 
O c Pn*I • reps 
A P M* 1 • VP* 

/Nil • IH 
AP'PNl • AC r PN 
Dl ?K?1-D r I’n? 

GG TC ? 

4 OryJ N p 

c ••••* CALL S'JsRruTlNtS Tu UP T * MiSTPFt AhO CR05S PLOTS* 

CALL T PI f T ( /flP *GK r 0A,CP5) 

CALL LPL^t | *AP ,UhFGA»CPS < 
r,n } 

CALL ’’l n T U.W< ,C#9V9) 
c Hl 


»♦* f N,r,y 1^1 

* QUALITY 



APPENDIX B - CONSIDERATIONS IN PILOT-INDUCED 
OSCILLATION SUPPRESSION FILTER DESIGN 


During the final flight of the shuttle approach and landing ♦ >, 

during the period when the pilot was experiencing a PIO, control uts 
were being made at a rate of 3 to 4 radians per second at peak amplitudes 
of ±10° of rotational hand controller travel. Since the inputs were 
sufficient to rate saturate the elevator, this experience was an excel' l 
example of the type of encounter a suppression filter should be designed 
to prevent. 

Figure 12 presents the rate saturation boundaries, the maximum 
surface travel limits, and the linear response of the elevator for 10° of 
controller travel, all as a function of frequency. The circle symbols 
represent nonlinear calculations of the actual rate limit response. The 
solid line indicates the total rate limiting that would occur in a 
sawtooth fashion. The dashed line represents the frequency and amplitude 
where barely perceptible rate limiting would occur. The PIO frequency, as 
pointed out, is between 3 and 4 radians per second. For example, if the 
PIOS filter is to avoid perceptible rate limiting, the forward loop gain 
would have to be reduced to approximately 38 percent of the linear value 
at 4 radians per second. 


Since the frequency is less than 10 radians per second, a rate 
- 1 (s + 20)^ 

5 (s + 10 r 


sensing element of 1 


was selected. This expression 


can be simplified to I s - ^^3) Figure 13 presents the 

41 (s + ior 

frequency response of this rate-sensing function ana illustrates the 
differentiation at low frequencies caused by the "free s" of the latter 
expression. As indicated, the gain is fairly linear wi n frequency from 
0 to 4 radians per second. Therefore, this function would be a reasonable 
estimator of frequency. The output of the rate-sensing function plus the 
control input is squared (fig. 2), then filtered with a first-order lag 
with a 3.3 second time constant. The smoothed output of the filter is 
then reverted back as the square root. The ratio of the two time- 

dependent parameters, ma Y thus be used as an estimate of the 

input frequency. EDEP. The F-l schedule (fig. 4) avoids division by zero 
and is further removed than necessary to give an unsuppressed step 
response for small controller inputs (inputs less than 4.5°). The F~2 
schedule (fig. 5) is used to limit the maximum amount of suppression. 
Selecting this value causes the output of EDEP to be determined as a 
function of frequency and amplitude (fig. 14). For a sinusoidal input, 
with a peak amplitude of 10° at the PIO frequency, the final value of 
EDEP would be approximately 0.3. In order to reduce the forward loop 
gain to 0.38 at 10° of controller input, the output of the pitch stick 
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shaping function would have to be reduced froa 6.98° to 2.62°. This will 
make the desired Minimal value 'or SKQ equal to -0.15. The schedule then 
is designed to intersect the ninimun value of EOEP at 0.3, as ;hown in 
figure 4(b). 
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Figure 2. PIOS filter logic and shuttle pitch stick shaping function. 
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Fiqure 4. PIOS filter schedules (fiq. 2). 


15 



C J 


o 1 2 3 4 t t 0 r -7 3 r r 6 “‘*7 

rime, sec Time, sec 

Figure 5. Output of the PIOS filter parameters with time for a sine wave Input. DEP - 10° 
sin cut (w = 3 rad/sec). 
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(a) DEP = 17° sin ut. 


Figure 6. Cross plot of the PICS filter characteristics through five 
complete sine wave cycles at a frequency of 3 rad/sec. 
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Figure 6. Continued. 
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(c) DEP = 5° sin ut. 


Figure 6. Ce icluded. 
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(a) 5° step input to DEP. 
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Figure o. Output to a step resonse in the pitch stick shaping logic 
with and without the P I OS filter. 
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Figure 9. Open-loop response of PIOS filter parameters with time. 
DEP, were taken from a shuttle approach and landing during which a 
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(a) Basic stick shaping. 


(b) With PIOS filter. 


Figure 10. Shuttle simulation using a closed loop attitude system with pure 
gain for a pilot model, K = 7.3 (see fig. 1). 
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Subsidence ratio 



(a) Variation of subsidence rat with pilot aain. 


Figure 11. Shuttle simulation usinq a closed-loop attitude system with pure 
gain for a pilot model (fig. 1). 
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(b) Variation of frequency with pilot gain. 
Figure 11. Concluded. 
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Figure 12. Elevator rate saturation boundaries and the predicted linear elevator response for 
the shuttle landing as a function of frequency for a pilot Input, DEP, of 10 degrees. 
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Figure 13. Amplitude ratio and phase angle, as a function of frequency. 
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